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Abstract 

We evaluate the neutral current quasi-elastic neutrino cross section within 
two nuclear models: the SuSA model, based on the superscaling behavior 
of electron scattering data, and the RMF model, based on relativistic mean 
field theory. We also estimate the ratio {vp — )■ vp) / {vN — > vN) and compare 
with the MiniBooNE experimental data, performing a fit of the parameters 
Ma and g^^ within the two models. Finally, we present our predictions for 
antineutrino scattering. 
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1. Introduction 

The study of neutral current mediated quasi-elastic (NCQE) neutrino- 
nucleus scattering in the GeV region is a powerful tool for hadronic and 
nuclear studies. We note that although in the tradition of neutrino exper- 
iments the term 'elastic', either neutral-current elastic or charged- current 
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elastic is used for neutrino scattering off free nucleons as well as on nucle- 
ons bound on nuclei, in this work we will refer to the latter case with the 
more precise denomination of quasi-elastic (QE). NCQE can be used, on one 
hand, to obtain information on the structure of the nucleon, in particular 
on its strange quark content, on the other it represents a probe of nuclear 
dynamics complementary to neutrino charged current quasi-elastic (CCQE) 
scattering and electron scattering. Several theoretical investigations have 
been devoted to the study of this reaction making use of different nuclear 

models QBBHHHBBi. 

The MiniBooNE experiment l3| has recently reported a high-statistics 
measurement of the NCQE cross section on mineral oil (CH2) and of the ra- 
tio {up — )■ up) / {uN — ^ uN) between single-proton and proton-|-neutron cross 
sections. In this letter we compare these measurements with the predictions 
of two relativistic nuclear models, the Super-Scaling- Approximation (SuSA) 
and the relativistic mean field (RMF) models, which have been previously 



applied to the CCQE process llll. Il2l |. A detailed description of the two 



models can be found in Refs. and [13[. Here we just recall their main 
ingredients: the SuSA approach is based on the assumption that the super- 
scaling function IJ] extracted from quasi-elastic electron scattering data can 
be implemented in the neutrino-nucleus cross section, the only differences 
between the two processes being related to the elementary reaction and not 
to the nuclear response; the RMF model provides a microscopic description 
of the process, where final-state interactions (FSI) are taken into account by 
using the same relativistic scalar and vector energy-independent potentials 
considered to describe the initial bound states. Both mode ls g ive an excellent 



representation of the experimental superscaling function [13|, in contrast to 
the relativistic Fermi gas (RFC), which fails to reproduce the electron scat- 
tering data. 

It has been shown in Ref. [l2| that, when applied to CCQE reactions, the 
RMF and SuSA models give similar results, although some difference arises: 
both models underestimate the MiniBooNE data 15|, but the RMF gives a 



smaller discrepancy. It has been suggested by various authors [16|, Il7|, Il8|, [19 
that the gap between theory and data can be filled by meson-exchange 
currents, multinucleon emission or particular treatments of final-state in- 
teractions. If one sticks to a simple nuclear description, such as the RFC 
model, presently used in neutrino interaction generators, the experimen- 
tal increase in the cross section can be obtained by introducing a nucleon 
axial mass Ma = 1-35 GeV, significantly larger than the standard value 
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Ma = 1.03 GeV [20||, which simulates the additional nuclear effects not con- 
sidered in the RFC 



2. Results and discussion 

Let us now consider the neutral current (NC) process. In order to compare 
with MiniBooNE data on CH2, we evaluate the following differential cross 
section per nucleon 



t;Cup,h{Q ) — — I- i=Cun,c{Q 



dQ^ 7 "^^'"'^ ' dQ^ 7— / ^Q2 

+ ^C^n,c{Q ) ^Q2 ' 

which results from three contributions: scattering on free protons, bound 
protons in Carbon and bound neutrons in Carbon, each of them weighted 
by an efficiency correction function Ci and averaged over the experimental 



neutrino flux 10| . Results corresponding to the two models mentioned above 
as well as the RFC are shown in Fig. [1 as functions of the "quasi-elastic" four- 
momentum transfer Qqe defined in 10| or of the outgoing nucleon kinetic 
energy T/v- The standard value Ma = 1.03 GeV has been taken for the 
axial mass, while the strange quark contribution to the axial form factor at 
= 0, Qa' (or equivalently As), has been set to zero. For the electric and 
magnetic strangeness the results of a recent global analysis of PV electron- 
proton asymmetry data [2H {ps = 0.59, /is = —0.02) have been used. Note 



however that the cross section is essentially independent of ps, fig |22 . 

We note that the SuSA cross section is smaller than the RFC one by about 
20% and the two curves have essentially the same slope in Q^. On the other 
hand the RMF result has a softer behavior, with a smaller slope. This 



is at variance with the CCQE case, for which, as shown in Ref. [12|, SuSA 
and RMF cross sections are very close to each other. This result indicates, 
as expected, that the NC data, for which the outgoing nucleon is detected, 
are more sensitive to the different treatment of final-state interactions than 
the MiniBooNE CC data, where the ejected nucleon is not observed. 

In Fig. |2]we illustrate the dependence of the cross section upon the axial 
mass Ma at strangeness Qa' = 0. We compare results with the standard 
axial mass to the ones obtained with the value of Ma that provides the best 
fit to the cross section within either SuSA or RMF models. We fit the axial 
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Figure 1: NCQE flux-avcragcd cross section computed using the RMF (solid blue), SUSA 
(dashed red) and RFG (dot-dashed green) models and compared with the MiniBooNE 
data 



mass performing a test using the true energy data from MiniBooNE 10 
(top panel on Fig. |2]) with the following definition 

Qgexp Qgtheo 



where CS^^^ is the experimental cross section in the i-bin, CS*'^'^° is the pre- 
dicted one and ACS^^'' is the error in CSfP. For g^^^ = 0, the l-a allowed 
regions of the axial mass for the two models are 

Ma = 1.34 ± 0.06 GeV for RMF (3) 
Ma = 1.42 ±0.06 GeV for SuSA, (4) 

corresponding to /DOF = 16.5/22 and x^/DOF = 4.7/22, respectively. 
These have to be compared with x^/DOF = 46.2/22 (RMF) and x^/DOF = 
45.3/22 (SuSA) for Ma = 1.03 GeV. 

In Fig. [2] the RMF and SuSA results are compared with the MiniBooNE 
data as functions of the true (top panel) and reconstructed (bottom panel) 
energies. Whenever a physical quantity is measured there are distortions to 
the original distribution in the observed quantity. Experimentalists correct 
the data distribution using unfolding techniques. There is an alternative 
method, which is to report them in the reconstructed nucleon energy, without 
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Figure 2: NCQE flux-avcragcd cross section computed within the RMF (solid blue lines) 
and SuSA (dashed red lines) models, compared with MiniBooNE data [10|] as a function 
of true energy on top panel and of the reconstructed energy on bottom panel, for different 
values of Ma (see text). 



applying the unfolding procedure (and corresponding errors). To produce 
the reconstructed energy results we used the folding procedure detailed in 



Appendix B of Ref. [23|. We observe that both models give a reasonably 
good representation of the data when the non-standard value of the axial 
mass is used. Moreover, we note that the SuSA cross section reproduces 
quite well the slope of experimental data, better than RMF one which has a 
smaller slope and falls slightly below the error bars for lowest data. It 
is however important to observe that none of the two models is expected to 
describe correctly the low-Q^ region, where collective effects play a dominant 
role. The values of the axial mass obtained with both models are compatible, 
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within 1-0", with the value Ma = 1-35 GeV employed by the MiniBooNE 
collaboration to fit their RFG model to the CCQE data. 



It has been known for some time (see, e.g., 2J, |25|) that the 



dependence of the NCQE neutrino-nucleon cross section is very mild. This 
results from a cancellation between the effect of g^^ on the proton and neu- 
tron contributions, which are affected differently by the axial strangeness: 
by changing g^f' from zero to a negative value the proton cross section gets 
enhanced while the neutron one is reduced, so that the net effect on the total 
cross section is very small. For that reason, the previous analysis performed 
to fit the axial mass is quite independent of the axial strangeness, which we 
just set to zero. Thus, once we have obtained the axial mass that fits the 
neutrino cross sections, we can look for a different observable that can be 
more sensitive to the axial strangeness content of the nucleon. Variations 
of the axial strangeness can have a large impact on the ratio between pro- 
ton and neutron cross sections. Furthermore, many systematic errors are 



canceled [26| in taking the ratio. 

The MiniBooNE experiment cannot measure the p/n ratio because the 
I'n — )• un reaction cannot be isolated. However, single-proton events can be 
isolated above the Cherenkov threshold, and so it was possible to construct 
two different samples: uN — uN (where N is either a proton or a neutron) 
with the standard NCQE cuts and a. up up NCQE proton-enriched sample 
for which two additional cuts were applied. The ratio {up — > up) / [uN uN) 
was reported in Ref. function of the reconstructed nucleon kinetic 

energy Tj-ec from 350 to 800 MeV. 

We now compare the predictions of our models with the experimental 



ratio, using the cross section folding procedure described in [23|. Following 
this procedure we convert our NCQE 'true energy' cross section into NCQE 
reconstructed energy distributions for the numerator and denominator sam- 
ples, separately, and finally we take the ratio. As expected, the ratio, unlike 
the cross section, is sensitive to axial strangeness. We now set the axial mass 
to the values (13|4|) previously found from the best fit of the cross sections at 
g^A = 0, and perform a fit to the axial strangeness parameter. The 1-a 
allowed regions turn out to be 

g^^^ = +0.04 ± 0.28 for RMF (5) 
^J) = -0.06 ± 0.31 for SuSA, (6) 

corresponding to /DOF = 33.6/29 and /DOF = 31.3/29, respectively. 
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In Fig. Owe present the ratio computed by using the above sets of param- 
eters as well as with the standard axial mass and no strangeness, as reference. 
From the comparison between the two curves in the top panel having stan- 
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Figure 3: Ratio {vp vp)/{vN — ?> vN) computed within RMF and SuSA models. Shad- 
owed areas represent the 1-cr region ahowed for (see text). The ratio computed with 

the hcsi-g\' is presented as weh as those obtained with the standard axial mass and no 
strangeness. Data from Ref. (loj . 

dard parameters Ma = 1-03 and g'^^ = (green and blue lines), it appears 
that the dependence upon the nuclear model is essentially canceled in the 
ratio, confirming that this is a good observable for determining the axial 
strangeness content of the nucleon. Within the error bars, the values of g^^ 
obtained are compatible with the ones of the previous analysis. Of course be- 
fore drawing definitive conclusions on the allowed value of g^^\ an extended 
analysis of the nuclear effects that are being effectively incorporated in the 



7 



S 0.01 . 
e 



0.001 




RMF;M^=1.03GeV;g^ =0 

— . SuSA; M^=1.03 GeV; g/'=0 

RMF;M^=1.34GeV;g/'=0.04 ; 

-- SuSA; M =1.42 GeV; g/"=-0.06: 



0.2 0.4 0.6 



0.8 1 



1.2 



1.4 



1.6 



Figure 4: NCQE antincutrino cross section computed using RMF and SuSA models for 
different values of and Ma- We employed the antincutrino flux prediction for Mini- 
BooNE given in Ref. 



increased value of Ma should be performed. However, it is worth mention- 
ing that the ratio shown in Fig. [3] shows little sensitivity to a possible np 
charge-exchange due to FSI. For instance a 20% of charge-exchange would 
not affect the results displayed in Fig. |3]by more than a few percent, for any 
reasonable value of g^^\ 

Before concluding we show our predictions for the NC antineutrino cross 
sections. In this case cross sections are slightly more sensitive to the axial 
strange content than neutrino ones, especially at high [iif. This is 

illustrated in Fig. |U where we display the flux-averaged NCQE antineutrino 
cross sections as a function of Qqe obtained within the RMF and SuSA 
models. We use the set of parameters ( 13]|4]|5]|6|) . favoured by our analysis of 
MiniBooNE NCQE neutrino cross section and ratio. As a reference, results 
for the standard axial mass and axial strangeness equal to zero are also 
presented. 

3. Conclusion 

Summarizing, we have presented predictions for MiniBooNE NCQE neu- 
trino cross sections with two nuclear models, SuSA and RMF. As it was the 
case of CCQE data, with the standard value of the axial mass, these models 
underpredict the cross section data. We have used the axial mass as an ef- 
fective parameter to incorporate nuclear effects not used in the models, such 
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as multi- nucleoli knockout. In doing this, we could restore agreement of the 
models to the data, with axial mass value consistent with the ones obtained 
in a similar fit to CCQE data by the MiniBooNE collaboration. The nu- 
clear models tuned this way can be employed to analyze NCQE cross section 
ratios as the ones obtained in the MiniBooNE experiment. We remind the 
reader that the present calculations are based on one-particle-one-hole as- 
sumptions and meson-exchange currents have not been considered. Provided 
that the models, as done here, are tuned to reproduce the total cross section 
data, our analysis shows that the ratio does indeed show a strong sensitivity 
to the axial strangeness content of the nucleon, while being highly model 
independent. This shows the relevance of having extended, good statistics 
data, eventually including antineutrino cross sections under similar condi- 
tions, which may help to disentangle the properties of the neutrino-nucleon 
and neutrino-nucleus interactions, of paramount importance for neutrino- 
oscillation experiments. 
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